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Abstract 

We study the large time behavior of non-negative solutions to the nonlinear diffu¬ 
sion equation with critical gradient absorption 

dtu — ApU + |Vm|^* =0 in (0, oo) x , 

for p G (2,oo) and g* := p — N/{N + 1). We show that the asymptotic profile of 
compactly supported solutions is given by a source-type self-similar solution of the p- 
Laplacian equation with suitable logarithmic time and space scales. In the process, we 
also get optimal decay rates for compactly supported solutions and optimal expansion 
rates for their supports that strongly improve previous results. 
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1 Introduction and results 


The goal of the present paper is to study the large time behavior of non-negative solutions 
to the following equation which combines a nonlinear diffusion and a gradient absorption 
term: 

dtu —/S.pU + \S/u\‘^* = 0 in(0, oo)xM^, (1.1) 

u{0) = uq in , (1.2) 

where the p-Laplacian ApU is given by ApU := div(|Vtt|P“^Vu) for p G (2,oo), and the 
critical exponent is g* := p — A^/(A^ + 1) > 1. We consider the Cauchy problem ()l.ip - (ll.2h 
for initial data uq satisfying 

Uq G 1T^’°°(M^), Uq >0, Uq ^ 0, suppU q is compact in M^. (1.3) 

The theory for the more general equation 

dtu — ApU + |Vu|'^ = 0, in (0, oo) x 

with general values of the exponents p and q, developed very quickly in the last years, after 
understanding better how the competition between the two effects—nonlinear diffusion in 
form of a p-Laplacian term, and first order absorption—works with respect to different 
ranges of these exponents. 

The semilinear problem, when p = 2, has been understood first, due to the possibility 
of using, at least partially, techniques related to the heat equation. It has been noticed 
that there exists two critical values of the exponent q, namely q = 1 and q = q* = 
(A^ + 2)/(A^ -|- 1), and the qualitative theory together with the large time behavior are now 
well understood after a series of works [DElEllHKinKISKle]. More precisely, for g > g* we 
have an asymptotic simplification with dominating diffusion, meaning that the absorption 
plays no role in the large time behavior, while for 1 < g < g*, there exists a special, unique 
self-similar solution of very singular type, which is the asymptotic profile for the evolution 
H]. The critical case corresponding to g = g* requires a different treatment and has been 
investigated in [H] using central manifold theory techniques. Finally, the case g = 1 is the 
subject of a well-known open problem (see [8] for the best estimates available), while for 
g G (0,1), finite time extinction occurs [U (TJ [15]. 

As heat equation techniques are not anymore available for p > 2, its study is more 
involved and it came later. In this range, the qualitative behavior is very different: on 
the one hand, the nonlinear slow diffusion implies a finite speed of propagation which 
entails that solutions emanating from compactly supported initial data stay compactly 
supported for all times. On the other hand, there exists a range of exponents g G (l,p — 1] 
where the dynamics is governed by the Hamilton-Jacobi part, a new fact that does not 
appear for p = 2 and g > 1; indeed, for g G (l,p — 1], the large time behavior is given 
by profiles in form of ” sandpiles” or ” regularized sandpiles”, reminding of the solutions to 
Hamilton-Jacobi equations, see [SUES]. For higher values of the exponent g, asymptotic 
simplihcation with dominating diffusion is expected for g > g* = p — N/(N -|- 1) and the 
existence of very singular self-similar solutions is known for g G (p — l,g*) [23], the lack 
of a uniqueness result preventing a complete understanding of the large time asymptotics 
in that case. In the present work we complete the panorama of the slow diffusion case by 
studying the case g = g*, where finer estimates for the solution than the general ones in 
[2] are needed. 
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We finally mention that, for p G (1,2), a qualitative theory is developed starting from 
the paper where a new critical exponent q = pj2 is discovered, and the qualitative 
theory for any g > 0 is established. In particular, there are again critical exponents q = q* 
and q = p/2, limiting ranges of parameters with different behaviors: diffusive for q > g*, 
algebraic decay for q € (g*,p/2), exponential decay for q = p/2, finite time extinction 
for q < p/2. Recent work by the authors helped to establish the existence of special 
solutions and the large time behavior for some of these cases, such as g G (g*,p/2), with 
very singular self-similar solutions, see [181120) . and q = p/2 where eternal solutions have 
been discovered m- The critical case g = g* is still open in this range. 


Main results. Coming back to Eq. our goal is to determine a profile that the 

solutions approach as t —>■ oo. As we are in a critical case which plays the role of an 
interface between purely diffusive behavior and mixed-type behavior, some logarithmic 
time scales are expected to appear. 

We introduce the following constant which will be important in the analysis 

' N{p-2)+p p{N + l)-2N ^ ^ 

We are now in a position to state our main asymptotic result. 


Theorem 1.1. Consider an initial condition uq satisfying q = q*=p-N/{N + l), 

and let u be the solution to the Cauchy problem (fTTTD-ra with initial condition uq. Then 


t^OO 

where 


lim 


u{t, x) — 


1 


(p-l)/(p-2) 


Ba, 


X 


t^(logt)(2-p)W+l)»7 


= 0 , 


and A* is uniquely determined and given by: 

I = - 


Bo = 

P 


„N-1 


dr 


p{p-2){N+l)ri/{p-l) 


V 






dr 


• (1-5) 


Let us remark that Ba is the well-known Barenblatt profile, the fundamental solution (and 
also the asymptotic profile) of the pure diffusion equation (without the gradient absorption 
term). The effect of absorption is seen in the fact that both the time-decay rate and the 
expansion of supports are changed with respect to the p-Laplacian equation. 

In order to prove Theorem 11.11 we have noticed that the precise time decay rate of the 
solutions was missing from the theory. Indeed, it follows from [2] that the L^-norm of 
the solutions to the Cauchy problem (HID-dOl) with initial condition uq satisfying m 
decays with a logarithmic rate as t —>■ oo; however, as we will prove, the rate in [2] is not 
optimal and it can be improved. Based on that, we find the exact logarithmic correction 
for the profile. Once optimal estimates are obtained, the final step will be an application 
of the stability technique of Galaktionov and Vazquez dais] together with some analysis 
concerning the variation of the L^-norm of the solution, which is needed to establish the 
uniqueness (and precise parameter A = A*) of the asymptotic profile. 


A by-product of our analysis is the following expansion property of the positivity set of 
non-negative compactly supported solutions to m- 
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Proposition 1.2. Consider an initial condition uq satisfying (1131), q = q* = p-N/{N+l), 
and let u he the solution to the Cauchy problem for the equation dHD with initial condition 
Uq. Introducing the positivity set 

V{t) := {x € : u{t,x) > 0} 

of u at time t>0, there are £»2 > £*1 > 0 such that 

B (o, ^it^(log c V{t) C B (o, ^2t''(log , t>l . 


Organisation of the paper. We begin by establishing our optimal time decay rate for the 
solutions to dnD, together with the rate of expansion of supports. As we already men¬ 
tioned, they strongly improve existing results and we consider them as the main theoretical 
novelty of the present work. This is done in Section [2l making use of a hne analysis of the 
supports and functional inequalities of Holder and Poincare type. The precise estimates 
are given in Proposition 12.11 and Corollary 12.31 These sharp estimates give in particular 
the correct time scales, thus allowing for a rescaling to reach a new equation with solutions 
being uniformly bounded in H The scaling step is performed in Sec¬ 

tion [31 where we also show that the Barenblatt prohle itself provides a suitable subsolution, 
which is fundamental in the sequel in avoiding the phenomenon of collapse to zero in the 
limit. Finally, in Section SI we complete the proof of Theorem 11.11 bv using the stability 
theorem for small perturbations of dynamical systems from nails], which we briefly recall 
in a more abstract framework in the Appendix. This approach is by now rather standard; 
nevertheless, we present it in some details, as the presence of the gradient absorption term 
leads to some further technical complications. 


2 Optimal L°°-bounds 

Let Uq be an initial condition satisfying (USD and denote the corresponding solution to the 
Cauchy problem (HU) by u. We recall that the existence and uniqueness of a non-negative 
(viscosity) solution u € i?C([0,oo) x to (11.11) which is also a weak solution follows 
from [21 Theorem 1.1], and moreover it satishes 

0 < u{t,x) < lluolloo, ||Vtt(t)||oo < IIVmoIIoo, 

for any t > 0. In addition, it follows from [U Theorem 1.1 & Theorem 1.2] and jH 

Theorem 1.6 & Corollary 1.7] that u{t) is compactly supported for any t > 0 so that the 

maximal radius of the (compact) support of u{t) and uq dehned by 

Q{t) := inf{i? > 0 : u{t,x) = 0 for |x| > R}, qq := q{0), 

are hnite for each t > 0. Furthermore, the following estimates for the support and the 

L^-norm are established in the above mentioned references; there is Ci > 0 such that 

Q{t) < C'i(l -|- t)'^ for any t >0 (2T) 

and 

||if(i)||i < Cl log(l -|- for any t > 1. (2.2) 
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However, as we shall see below, the estimates m and (12.2p are far from being optimal 
and can be strongly improved. We finally recall that, due to [H Lemma 4.1], the solution 
u to dni) cannot vanish identically in finite time, that is, 

u{t) ^0 for t> 0 . 

Proposition 2.1. There exists a positive constant C 2 , depending only onp, q, N and uq, 
such that for any t > 0, 

Qit) < C2(l+t)'^(log(l + t))-’'(P-2)(^+i), (2.3) 


and 

\\u{t)h<C 2 {log{l + t))-^^+^\ (2.4) 

Notice that this is a real improvement in the second estimate (12.4p with respect to (|2.2p . 
since 

1 _ A^ + 1 ^ A^ + 1 _ Yll 

g* - 1 ^ p(iV + 1) - (2iV + 1) 2{N + 1) - (2iV + 1) ^ 

for p > 2. 

Proof. Since uq is non-negative, continuous and compactly supported in M^, there exists a 
non-negative, continuous and radially symmetric, radially non-increasing function Uq with 
compact support, such that 0 < uo(x) < Uq{x) for any x € M^. Let U be the corresponding 
solution to (11.11) with initial data Uq] it follows from [SJ Theorem 1.2] that the function 
X U{t, x) is also radially symmetric, radially non-increasing and compactly supported, 
for any t > 0, and we deduce from the comparison principle that 0 < u{t,x) < U{t,x), for 
any {t, x) € (0, 00 ) x Let, for t > 0, 


a(t) := sup {ja;| : U{t,x) > 0} 


be the radius of the support of U. We then have g^t) < cr{t). Furthermore, we infer from 
[21 Theorem 1.2 &: Proposition 1.4] that 

|Vt/(P-2)/(p-i)(t,a;)| < C3|lt/(s)|l^-2)/P(P-i)(t - s)-^/^ for 0 < s < t (2.5) 

and 

||t/(i)||oo < Cs\\U{s)\\Y>{t - s)-^\ for 0 < s < L (2.6) 

It next follows from (|l.ll) that, for any non-negative function y € (^^([O, 00 )), we have 


dt 


{l^l>y(t)} 


U{t, x) dx = 


'{\^\>y(t)} 

[ 


dtU{t, x) dx — y'{t) 


K\A=y{t)} 


U{t, x) dx 


< I div(jVC/|^ '^VU){t,x) dx — y'{t) / U{t,x)dx 


'{\A=yd)} 


<- I \VU{t,x)\P-^VU{t,x) ■ ^dx 

J{\x\=y{t)} fI 

— y'{t) / U{t,x)dx 

d{\x\=y{t)} 

> - 1 


< 


'{\x\=y{t)} - 


p — 2 


\vu^p 


-2)/(p-1)f 


{t,x)\P -y'{t) 


U{t, x) dx. 
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Fix now to > 0- For t > to we deduce from the estimates (12.51) (with s = (t + to)/2) and 
(j2.6p (with ((t + to)/2, to) instead of (t,s)) that 


p — 1 
p — 2 


|VC/(P-2)/(P-F(t,x)|P-^ < —lc7P-i u ( 

p — 2 y 


t + to 


P-2)Ip 




< 


(2Cf)(P"F/p^ 


- 1 


c 3 \\u{tom 


PV 


P-2 

<vC4uito)\\^f-^'^^{t-tor 


t - to 


-Npl (P-2)/P 


- +„i-(P-l)/P 


{t - to) 


-1 


Choosing 

y{t) := cj(to) + C 4 \\U{to)\\‘f~‘^'^'^{t - to)^, 
the above inequality reads 

^|VC/(P-2)/(P-F(t,x)|P-i<y'(t), 

p — 2 

for any t > to ^ from which we deduce that 

d f 

— U{t,x)dx<0 for t > to- 

at J{\x\>y{t)} 

Since y{to) = a (to), we end up with 


/ U{t,x)dx < / U{to,x)dx = 0, 

J{\x\>a{to)} 

for any t > to- Owing to the non-negativity of U, this is only possible if U{t,x) = 0 for 
|ic| > y{t), which means 

cr{t) < y{t) = o-(to) -bC'4||[/(to)||i^"^^''(t - to)^, (2.7) 

for t > to- In order to proceed further, we need the following Poincare inequality: 

Lemma 2.2. Given y G [1, oo) and R > 0 there exists a constant K depending only in N 
and y such that 

for all w €Wo’^{B{0,R)) . 

Proof of Lemma \2.2[ This follows by combining Holder and Poincare inequalities. More 
precisely, setting y' := y/{y — 1), one has 

where we have used that the Poincare constant is of order 0{R). □ 
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In order to continue the proof of ProDosition l2.ll we first note that integrating (jl.ip over 
gives 


A 

dt 


rwiii + iivc/|i^: = o. 


Since U(t) is supported in B{0,a{t)) for each t > 0, it belongs to cr(t))) and 

we apply Lemma 12.21 to obtain 


, 1 \mmr .. 

for any t > 0. Now, fix T > 1 and introduce the notation 

S(r) := sup {f^(logt)^cr(t)} , A:=r]{p-2){N + 1). 

te[i,r] 


( 2 . 8 ) 


We infer from ()2.8p that, for any t G [1,P], 


whence 


dt" " Kt {t-^{\ogt)Aa{t)Yh 


^ li-w.M, (logt)^^ iirr^.tMI<?* 


< 0 , 


+ Kt{T.{T)Yh 


\\U{t)\\f<0, te[l,T]. 


Integrating the above inequality over (l,t), t G (1,T), we hnd 

||C/(t)||i < {K{N + l))i/(''*-i)S(r)i/'?(5*-^)(logt)-(^+^\ t G (1,T]. (2.9) 

Next, let m > 1 be an integer to be determined later and consider t G (1,P]- 

• Either t < 2™ and it follows from (12.71) with to = 1 that 

t-'?(logt)^a(t) < (log2"*)^ [a(l) +C 4 ||t/( 1 )||S'’""^^ 2 ™^] < C(m) . 

• Or 2”^ < t < T (if 2"* < T), and we infer from ()2.7p (with to = t/2) and ()2.9p that 


t-^(logt)^u(t) < t-''(logt)^a ( - ) + C 4 t-''(logt)^ 


U 


< 2 " 


/ logt Y 

Vlog(i/2)y' 


log hr o 


{p-2)v / 

2 , 


+ C' 42 -’^(log t)^ j:{TY^viq*-l) I log 

Since t > 2™, we obtain that 


_(Ar+l)1 (p-2)»7 


logt < 


m 


log ( ^ ) , 


m — 1 

whence, plugging this estimate in the previous inequality, we get 

A 


t '^(logt)^cj(t) < 


m 


m — 1 


2-v 


s(T) + C4S(r)ip-2)/(-?*-i) 
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Combining the previous two estimates and taking the supremum over t € (1, T], we obtain 

S(T) < f —!^) 2-^S(T) + + C(m). 

\m — 1/ 


We now hx m large enough such that 

■ — 


m 


m — 1 


2-^ < 1 


Taking into account that p — 2 < g* — 1, we deduce from Young’s inequality that 

S(r) < 6m^{T) + i^S(r) + C(m) , 

which readily implies that S(T) < C{m) for each T > 1, the constant C{m) being inde¬ 
pendent of T. We have thus proved that 

a{t) < Cf^ilog ^ i 


while the uniform bound for S(r) together with (12.91) imply 

Recalling that u{t,x) < U{t,x) for any {t,x) G (0, oo) x and that g{t) < a{t), we 
obtain the expected estimates ([ 23 ]) and ([ 23 ]). □ 


Corollary 2.3. There exists a constant Cs > 0 depending only on p, q, N and uq, such 
that for any t > 0, we have 

Mt)\\oo<C,{l + t)-^^log{l + t))-P^^^+^^ ( 2 . 10 ) 

and 

||Vu(t)|U < C5(l + t)-('^+')"(log(l + (2.11) 


Proof. We combine the estimates in [21 Proposition 1.4] with the previous estimates of 
Proposition 12.II We thereby obtain 


k(t)||oo < C 




prj 


^—Nrje^Nri 


/ .\-{N+l)pr) \ -(Ar-|-l)p ?7 

< Ct-^^ log (^1 + - j < Ct-^^ (^- log(l + t)j 

< C(1 + t)-^^(log(l + t))-(^+i)P^ 


and 


||Vu(t)||oo <c 


u 



2r, 

^-n{N+i) 

1 


< C't-^(^+^)(log(l + 


which completes the proof. □ 


As we shall see in the next section, these estimates are optimal and are building blocks 
in identifying the large time behavior of compactly supported solutions to dni) and thus 
proving Theorem ll.il 













3 Scaling variables 


Let uq be an initial condition satisfying m) and denote the corresponding solution to 
the Cauchy problem ([nD-dOD by u. According to the estimates derived in Section [21 we 
introduce the following new variables (s,y) and function w: 



s = log(e + t), y = x{e + t) '^log{e + t)^^ 2 )(Ar+i)? 7 ^ 

(3.1) 

and 

u{t, x) = (e + log(e + t)~P'^^^^^'>w{s, y). 

(3.2) 

By (11.ip the rescaled function w solves 



dgW — Cw = 0 in (1, 00 ) x , 

(3.3a) 


rc(l) = Uq in , 

(3.3b) 

where 

Cz \= rjy ■ Wz + pNz + ApZ 



1 

(3.3c) 


-Vz — pri{N + l)z + {p — 2)7]{N + l)y ■ Vz] . 

s 

For further 

use. we introduce the autonomous counternart of (l3.3aD which is 



dgV — Lv = 0 in (1, 00 ) x , 

(3.4a) 

with 

Lz -.= gy - Vz + i]Nz + ApZ. 

(3.4b) 


The boundedness of w readily follows from the previous section. 


Lemma 3.1. There is a positive constant Cq depending only on p, q, N, and uq such that 

ll'«^(•s)l|l + ||'i«(s)||oo + ||V'u;(s)||oo < Ce, 

for any s > 1. Moreover, the support of w{s) is localized: there exists Rq > 0 such that 
supp('u;(s)) C B{0,Ro) for any s G (l,oo). 

Proof. The bounds for the lT^’°°-norm of w{s) are immediate consequences of estimates 
(|2.10D and (|2.11|1 . taking into account the definition of w in (13.2p . The estimate for the 
L^-norm follows from (|2.4p by a change of variable as below: 

\\w{s)\\i = {e + t)^^log{e + t)P^^^+^^ [ u(t,y{e + t)nog{e + t)-^P-'^^^^+^^^) dy 

Jrn \ J 

= log(e + t)^+^l|'u(t)||i < C 2 . 

Finally, the assertion about the localization of the support follows from estimate (12.3p and 
the definition of the new variable y in (|3.ip (with Rq = 02 ). □ 

Lemma 13.11 provides a fine upper bound for w but its optimality can only be guaranteed 
by a lower bound of the same order. In addition, such a lower bound would prevent the 
possibility of collapsing to the trivial solution in the limit s —>■ 00 . This is a consequence 
of the following result. 
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Lemma 3.2. There exists Agub > 0 sufficiently small such that the Barenblatt profile 

BA{y) = (A- , Bo = 

V / + p 

is a subsolution to Eq. ()3.3ap for A £ (0,^s„b). 

Proof. A simple computation shows that LBa = 0 for any A > 0, where L is the au¬ 
tonomous operator defined in 13.4bp . Moreover, 

VBA{y) ■ y = (a - Bo\y\P/^P-^A'^^"~"^ 

p — 2 \ / + 

hence 


-sCBa = + r}{p -2){N + l)y ■ VBa - r]{N + 1 )pBa 

- ffiN + l)pBo (a - Bo\y\P/^P-^^y^^''~^^ 

- ffiN+i)p (a - 

- ffiN + l)pA (a - Bo\y\P/^P-^^y^^~^^ 

= y(^A-Bo\y\P/^P-^'>yy"~"^ 

^(g*-p+i)/(p-i) (A-Bo\y\P/^P-^^y'^*~^^^^^~^^ |y|9*/(p-i) _ (N + l)pA 


Either |y| > {A/Bo)^p ^'>/p and —sCBA{y) = 0. Or \y\ < (A/iJo)*'^ and, since A 
Bo\y\P/'^P~^'^ < A, we find 


— sCBa 

<y{A-Bo\y\P/^P-^^) 

<r?A 
where 


l/(p-2) 

+ 

l/(p-2) 


^(<?.-p-Hl)/(p-l)^(g,-l)/(p-2) ^ - (A^ + l)pA 

^(q.-p+l)/(p-l)^-9*/p^0 - (AT + l)p 


Q g* - 1 (A^ + l)(p - 1)((A^ + l)p - 2N) ^ ^ 

p — 2 p p{p — 2) 

for p > 2, taking into account that q* = p — N/{N + 1). Consequently, there exists 


Asub 


{N + l)pBQ^^r] 


ni/e 


> 0 , 


such that —sCBa < 0 in (l,oo) x for A £ (0, A^^;,), ending the proof. □ 


10 



















A first consequence of Lemma 1,3.21 is the optimality of the temporal decay estimates 
established in Lemma l3.11 

Proposition 3.3. There are Aq £ (0, Agub) cind to > 0 such that 


n(i)||oo > Ao(e + t)-'^’'log(e + , 

t > 0 , 

(3.5) 

uito,x) > Bao{x) , X £ . 


(3.6) 


Proof. Since uq ^ 0 and is continuous by (ll.3p . there are xq £ , qq > 0, and Eq > 0 

such that 

uo{x) >eo, X e B{xo, qq) . (3.7) 

Introducing Uo{x) := Uo{x + xq), x £ and denoting the corresponding solutions to 

(jl.ll) and (I3.3a|) by u and w, respectively, we observe that the invariance of (II.1|) with 
respect to translations entails that u{t,x) = u(t,x + xq), while we infer from (13.711 that 

uo{x) > Ba{x) , X £ , (3.8) 

provided 

A < 4p-2)/(p-i) and A < BoqI'^^-^'^ . (3.9) 

Choosing A £ (0,As„fe) satisfying (13.91) . we deduce from Lemma [221 (13.81) . and the com¬ 
parison principle that 

w{s,y) > BAiy) , (s,y) £ (l,cx)) X . 


Coming back to u, we realize that 


(e + t)^^ log(e + t)P(^+b^ u{t, x) > Ba 



Xo) log(e-|-t)(^ 2)(Ar+l)?7 
(e + 1 )^ 


(3.10) 


for {t,x) £ (0,oo) X M^. A first consequence of (I3.10p with x = xq is the lower bound 
(I3.,^p . We also infer from (I3.10p with x = 0 that 


(e + t)^^ log(e + u{t, 0) > Ba 


( (xo) log(e-|- t)*^^ 2)(Af-|-l)r) 

(e + t)^ 


and the right-hand side of the above inequality is positive provided t is large enough. 
Therefore there is to > 0 such that u(to, 0) > 0, and we argue as in the proof of (13.8p to 
complete the proof of (j3.6l) . possibly taking a lower value of A if necessary. □ 


4 Convergence. Proof of Theorem 11.1! 

Thanks to the outcome of Sections [2] and [3| we are in a position to prove Theorem 11.11 
To this end, we follow the lines of the analysis developed by Galaktionov &: Vazquez in 
nans], the central tool being a stability theorem which is recalled in Section for the 
reader’s convenience. 


11 












We fix an initial condition uq satisfying m and denote the corresponding solution to 
(frT])-(fL2D by u. By Proposition 13.31 there is to > 0 and G (Oj^suh) such that 

u{tQ,x) > Baq{x) , X G . (4.1) 

We then dehne w by (ISID-B with u{- + to) instead of u, that is, 

u(t + to, x) = (e + t)"^*^ log(e + y) , 

the variables (s, y) being still given by (13.11) . We infer from (14.ip . Lemma l3. 11 Lemma 13.21 
and the comparison principle that, for all s > 1, 

||rc(s)||i + ||u;(s)||oo + ||Vt(;||oo < Q , (4.2) 

and 

w{s,y) > BAoiy) , y G , and w{s,y) = 0, \y\ > Rq ■ (4.3) 

We define the set 

X ■={z£ L^(M^) n BC{R^) : z{y) > BA^iv) , y G and ||z||i < Cq] , 

which is a complete metric space for the distance induced by the L°°-norm, the parameters 
Co and ^o being given in (|4.2I) and ()4.3p . respectively. We also set 

S := {rc} , 

and deduce from (13.ip . (I4.2p . (|4.3I) . and the properties of u that w G CPPO, oo); X). We now 
check that the set S enjoys the three properties required to apply the stability 

result from [121 [T3] recalled in Theorem lA.ll below. In our setting, the non-autonomous 
operator is the operator £ defined in (|3.3cl) and its autonomous counterpart L is defined 
in (j3.4bp . the associated evolution equations being (I3.3ap and (I3.4ap . respectively. 

Clearly, tc is a solution to (I3.3ap and it readily follows from the W^’°°-bound (14.21) . the 
uniform localization of the support (14.3p . and the Arzela-Ascoli theorem that {tc(s)}s>o is 
compact in so that (HI) is satisfied. 

We next infer from the same properties (14.2p and (14.3p that, for {s,y) G (0,oo) x , 
\£w{s,y) - Lw{s,y)\ < j [||Vu;(s)||^ + ||u;(s)||oo + |y|||Vu;(s)||oo] 

< j [Ct +C6 + RoCe] . 

Therefore, 

||£u;(s) - Lu;(s)||^ < ^ , s>0, 

from which (H2) follows. 

Finally, we fix Ai G (Ao,oo) such that Ai > BqRq^^ and ||i?Ail|i > Cq and set 

n:={BA : ^ G [Ao,Ai]} . 

Clearly is a non-empty and compact subset of X. Since (H3) concerns only the au¬ 
tonomous operator (j3.4aD . which also arises from the standard p-Laplacian equation via 
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a self-similar change of variable, the uniform stability of n holds true as a by-product of 
classical results for the p-Laplacian equation (see m Section 6] and also [131 Section 4.6]). 

We have thus checked the validity of (H1)-(H3) and may thus apply the stability the¬ 
orem (Theorem lA.lD to conclude that the cj-limit set (for the topology of the uniform 
convergence) of w is included in 11, that is, 

, , _ J u) G X : there is a sequence (sj)j>i ,Sj^oo, 

^ ^ ' I such that limj_^oo ll^(sj) “ halloo = 0 



Mass analysis. Uniqueness of the limit. It remains to show that the asymptotic 
profile is in fact a uniquely determined Barenblatt profile with the parameter A* as in 
Theorem ll.il To this end, we perform a mass analysis along the lines of a similar argument 
in [12]. We first observe that classical properties of the Barenblatt profiles ensure that, 
given s > 0, there is a unique A(s) > 0 such that 


©(s) := ||n;(s)||i = . 


(4.5) 


More precisely, 

0(5) = , (4.6) 


where 

h 

and ujn denotes the 
Since 


:= {N + 1) (l - dr (4.7) 

volume of the unit ball of M'^, see [241 Section 11.4.1] for instance. 


II^Ailli >Cq> 0(s) > 1|Bao1|i , s > 0 , 

by ()4.2I1 . (14.31) . and the choice of Ai, we deduce from the monotonicity of A i—)• UByijji and 
(j4.5l) that 

^0 < A{s) < , s > 0 . (4.8) 

In addition, integrating (I3.3a|) with respect to space shows that, for s > 0, 


^ ^ith G{z) := {N + l)||zlli - llVzji^: . (4.9) 

©wing to (|4.9p . (|4.6I) . and the regularity properties of w, we realize that A G C'([0, oo)) fl 
C^((0,oo)). 

We next claim that 

lim ll'u;(s) - = 0 . (4.10) 

s—^oo ^ ' 

Indeed, assume for contradiction that there are an increasing sequence {sj)j>i of positive 
real numbers, Sj —)■ cx), and s > 0 such that 


||u;(sj) - ^A(s,-)l|oo > e , j>l- (4.11) 

0n the one hand, we infer from (USD, (031), (03D, and the Arzela-Ascoli theorem that 
there are a subsequence of {sj)j>i (not relabeled) and A G [Aq, Ai] such that 


lim lltc(sj) — B^Woo = 0 . 
j-foo 
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On the other hand, it follows from ()4.8p that, after possibly extracting a further susbe- 
quence, we may assume that there is A^o G [Ao,Ai] such that —>■ Aoo as j —)■ oo. 

This readily implies that 

lim ||fiA(s,) - ^^oolloo = 0 . 

We may then let j —)■ oo in (|4.11l) to conclude that 

II^A-^A„elloo >e . (4.12) 


Now, 


= lim 0(sj) = lim 

j—>-oo J—>-oo 


B 


A{sj) 


= \\Ba. 


so that A = Aoo , which contradicts (I4.12p . We have thus proved (I4.10p . 

We next infer from ()4.2I) and (jd-.^p that, introducing 

/(s, y) := yy • Vw{s, y) + 7]Nw{s, y) 

- - [|Vrt;(s,y)|'^* - pr]{N + l)w{s,y) + {p - 2)7]{N + l)y ■ Vw{s,', 


for {s,y) G (0, oo) x Eq. (j3.3al) reads 


dgW — ApW = f in (1, oo) x i?(0, i?o + 1) 

with w G oo; i?o + 1))) and / G L°°((l,oo) x B{0,Ro + 1)). We then 

infer from m Theorem 1.1] that there are Cg > 0 and v G (0,1) such that 

|Vr(;(si,yi) - Vw{s2,y2)\ < Cg (^\yi - y2\'' + |si - S2\''/‘^'^ 


for all S 2 > Si > 2 and (yi, 2 / 2 ) G B{Q, Rq) x 5(0, Rq). Combining the above property with 
we deduce that {Vt(;(s)}s>i is bounded and equicontinuous in C{B{0, Ro)-,'R^) and 
thus compact in that space by the Arzela-Ascoli theorem. Recalling ()4.3p and (|4.10p we 
conclude that 


(I|w^(s) - 5A(s)||r + ||Vu;(s) - = 0 , rG[l,oo]. (4.13) 

An immediate consequence of (I4.13P is that 

^li^ |G(u;(s)) - G(5^(^))| = 0 , (4.14) 

while the explicit formula for gives 

= g{Ais)) (4.15) 


with 


9 {a) := , a > 0 , 

where R is defined in (|4.7I) . and 

I 2 := J ^((p-i)(^-i)+9)/(p-i) , 
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We observe that g vanishes only once in (0, cx)), for a = 
reads 



p(p-2)(Af+l)r?/(p-l) 


which is defined in (II. and 


with the notation introduced in this section. In fact, 

g{a) <0 for a > A^, and g{a) >0 for a G (OjA*) . (4-16) 


Thanks to (14.4p . (14.9p . (I4.13p . (I4.14K . (I4.15K . and (I4.16p . we are in a position to proceed as 
in the proof of m Proposition 5.2] to establish that uj{'w) = {Ba^}- Undoing the rescaling 
(I3T1)-(I321), this property readily gives Theorem 11.11 

We finally note that Proposition 11.21 follows at once from (14.3p after translating these 
properties in terms of u. 


A The stability theorem 

We briefly recall here for the reader’s convenience the stability theorem introduced by 
Galaktionov and Vazquez in [la US] and used in Section 01 As a general framework, 
consider a non-autonomous evolution equation 

ds^ = Ld , (A.l) 

that can be seen as a small perturbation of an autonomous evolution equation with good 
asymptotic properties 

ds^ = , (A.2) 

in the sense described below. We consider a set S of solutions i? € CdO, oo);X) to (lA.lK 
with values in a complete metric space {X,d). We assume that: 

(HI) For each G <S, the orbit {i?(t)}t>o is relatively compact in X. Moreover, if we let 

id^it) := '&{t + r) , t >0 , r > 0, 

then {i?’^}r>o is relatively compact in (7(10,T];V) for any T > 0. 

(H2) Let I? G 5 for which there is a sequence of positive times {tk)k>i, ik oo, such that 
+ ik) —^ 0 in C'([0,T]; V) as k ^ oo for any T > 0. Then 0 is a solution to 

(M). 

(H3) Define the w-limit set D of (IA.2I) in X as the set of / G A enjoying the following 
property: there are a solution <I> G (7(10,00); A) to ()A.2p and a sequence of positive 
times {tk)k>i such that —>■ oo and — > f in A. Then D is non-empty, 

compact and uniformly stable, that is: for any e > 0, there exists <5 > 0 such that if 
4* is a solution to (|A.2p with (i($(0),U) < 5, then (i(4>(t),U) < e for any t > 0. 

The stability theorem (also known as the S-theorem) then reads: 

Theorem A.l. If (H1)-(H3) above are satisfied, then the uj-limit set of any solution 
D (z S is contained in U. 

For a detailed proof we refer the reader to nans]. 
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